Pure Spin Currents and associated electrical voltage 
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We present a generalize Landauer-Biittiker transport theory for multi-terminal spin transport in 
presence of spin-orbit interaction. It is pointed out that the presence of spin-orbit interaction results 
in equilibrium spin currents, since in presence of spin-orbit interaction spin is not a conserved quan- 
titative. Further we illustrate the theory by applying it to a three terminal Y-shaped conductor. It 
is shown that when one of the terminal is a potential probe, there exist noneqmhbrium pure spin 
currents without an accompanying charge current. It is shown that this pure spin currents causes 
a voltage drop which can be measured if the potential probe is ferromagnetic. 
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Producing and measuring spin currents is a major goal 
of spintronics. The standard way is to inject spin currents 
from a Ferromagnet into a semiconductor in a two termi- 
nal geometry [1]. However this has a drawback, due to 
conductivity mismatch, the polarization of injected cur- 
rent is rather small and it always has an accompanying 
charge current [2], Also for any spintronics operation 
spin orbit interaction plays an important role, for e.g., in 
Datta-Das spin-transistor [3]. 

In light of these development it would be interesting 
and highly desirable if one can produce spin currents in- 
trinsically. One such possibility is provided by intrinsics 
spin-orbit interaction. Presence of impurity atom or de- 
fects gives rise to spin-orbit interaction of the form [4] 
[5], [6], 



H so = \{VU{r) x k) • tr 



(1) 



where a is a vector of Pauli matrices ,U(r) is potential 
due to defects or impurity atoms and k is the momentum 
wave vector of electrons and A is spin-orbit interaction 
strength. For strictly two dimensional case for which the 
potential U(r) depends on x and y coordinates only the 
Hamiltonian commutes with o~ z , hence z component of 
the spin is good quantum number. As is well know that 
this kind of spin-orbit interaction has a polarizing effect 
on particle scattering [7], i.e, when an unpolarized beam 
is scattered it gets polarized perpendicular to the plane 
of scattering. Further scattering of this polarized beam 
causes asymmetry in scattering processes, i.e. electrons 
with one particular spin direction, e.g., spin-up electrons 
have a larger probability to be scattered to the right com- 
pared to spin-down electrons [6], [7]. This property of 
spin-orbit scattering gives rise to novel effects like spin 
hall effect [4]. 

Here in this article we show that the the above dis- 
cussed property of spin-orbit scattering can be used to 
generate and measure spin-currents [5]. Consider a three 
terminal, two dimensional Y shaped conductor shown in 
Fig. 1 . The plane of conductor is xy. Since the conductor 
is two dimensional which fixes the scattering plane, the 
scattered electrons will be polarized along z axis (per- 
pendicular to the scattering plane). However the polar- 



ization for the two branches of Y junction will be oppo- 
site [4]. Hence a three terminal structure would create 
spin currents from an unpolarized current in presence 
of spin-orbit interaction [6] , [8] . Moreover a three termi- 
nal device provides an important possibility of generating 
nonequilibrium pure spin currents without an accompa- 
nying charge current. This is the case when one of the 
terminal acts as a voltage probe. For e.g., say the termi- 
nal 3 is a voltage probe as shown in Fig.l, i.e., the voltage 
V3 at third probe is adjusted such that the total charge 
current flowing in terminal 3 is zero, i.e. 1^=0 [9] , [10]. 
Physically it implies that the charge current flowing in 
(which is polarized as argued above) is same as charge 
current flowing out. However the polarization of charge 
current flowing out need not to be same as the polariza- 
tion of charge current flowing in(see Fig. 1). Hence there 
will be a net spin current flowing without the accompa- 
nying charge current. This is a pure nonequilibrium spin 
current. 

We support our prediction by generalize Landauer- 
Biittiker charge transport for multi-terminal spin trans- 
port. We provide unambiguous definition of spin cur- 
rents. Using this theory we discusses the possibility of 
generating and detecting nonequilibrium pure spin cur- 
rents and point out the existence of equilibrium spin cur- 
rent. 

The possibility of injecting pure spin currents were first 
discussed in Ref. [12] for a three terminal device where 
two of the terminals were ferromagnetic. Also direct opti- 
cal injection of pure spin currents in GaAs/AlGaAs quan- 
tum wells was demonstrated in Ref. [13]. We would stress 
that in our case spin current is not injected rather gener- 
ated intrinsically due to the spin-orbit interaction with- 
out any magnetic element in the system, which is not the 
case in Ref. [12]. Thus we avoid the problem of spin injec- 
tion altogether. Further, since the effect discussed relies 
on the general scattering properties due to the spin-orbit 
interaction. Hence it will be observable with any kind of 
spin-orbit interaction, e.g.. Rashba spin-orbit interaction 
[11]. 
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FIG. 1. Y-shaped three terminal junction with applied 
voltages Vi, V2 and V3 as depicted. Terminal third (labeled 
3) is a voltage probe (non-magnetic or ferromagnetic) which 
draws no charge current. However the polarization incoming 
and outgoing electrons are opposite to each other, causing a 
pure spin current 

We first briefly outline the spin transport theory for 
multi-terminal devices. Let us consider the two dimen- 
sional Y-shaped structure shown in Fig.l. The plane of 
structure is xy and a perpendicular to it defines the coor- 
dinate system. Let us choose the spin quantization axis 
to be along u, pointing along (8, <f>) where 6 and are 
usual spherical angles ( In other words we choose the 
spin basis to be eigen states of operator <x ■ u ) This is 
essentially since a charge current flowing along a spatially 
direction can be polarized along a direction which need 
not coincide with the direction of flow of charge current. 
Also in presence of spin-orbit interaction the rotational 
invariance in spin space is lost [14] , hence any theory for 
spin transport should take this fact into account. With 
this definition we can generalize Landauer-Biittiker the- 
ory for spin transport. Let V m be potential at terminal 
m measured from the minima of lowest band, where m 
can take values 1, 2 and 3 corresponding to the three 
terminals of Fig. 1. is spin resolved transmission 

probability of electrons incident in lead m in spin channel 
a to be transmitted into lead n in spin channel a. We 
point out that a and a need not to be same in presence 
of spin-orbit interaction, since SO interaction will make 
spin flip transmission probability non zero. The spin cur- 
rent a flowing into terminal m is, (here a can be either f 
or I ) 



r = — 



\ ± nm v m - L mn v n) 



(2) 



where a and a are indices labeling the two spin eigen- 
states for a chosen quantization axis. In writing above 
equations we have made an assumption that the spin re- 
solved transmission coefficient arc energy independent. 
A generalization of the above equation when the spin 



resolved transmission coefficient are energy dependent 
is straight forward. It amounts to replacing T"^ by 

Since SO interaction preserves time reversal symme- 
try, which lead to the following constrains on the spin- 
resolved transmission coefficient, 



rpa a 
- 1 nm 



(3) 



Using cq. 1 we can immediately write down the net 
charge and spin current flowing through terminal m, 



79 = r 
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where is charge current and 1^ is spin current. We 
stress that eq. (5) correctly determines spin current gen- 
erated by presence of spin orbit interaction. Since in 
the absence of spin-orbit interaction and any magnetic 
element in the device , spin resolved transmission coef- 
ficient obey a further rotational symmetry in spin space 
i.e. T"^ = T~^~ a , which implies that spin currents are 
identically zero for all terminals, i.e., 1^ = 0. 

Equilibrium spin current : To discusses equilibrium 
spin currents let us consider the case when all the po- 
tential are equal i.e., V m — Vq V m. In this situation 
charge current flowing in any terminal should be zero 
(J9j = 0) which leads to the following sum rule (from eq. 
(4)) 



^ ^ m — ^ ^ T^rr 



(6) 



where T nm — J2 a a ^nm * s total transmission probabil- 
ity(summed over all spin channels) from terminal m to 
n. This sum rule is robust and should be satisfied irre- 
spective of the detailed physics [15]. This is a well known 
gauge invariance condition. Charge conservation implies 
J2 m -^m = which follows from symmetry of spin resolved 
transmission coefficient ,eq. (3), and the gauge invariance 
condition, eq.(6). So in equilibrium there are no charge 
currents flowing. However this is not the case for the spin 
currents. This point can be appreciated if we look closely 
at the equation (5) for spin current. Since in general the 
transmission coefficient, T"^ ^ Enm* > which occurs in 
eq. (5). Hence even when all the potential are equal , 
the spin current given by eq.(5) is non zero. This is equi- 
librium spin current. Notice that this is consistent with 
time reversal invariance (eq. (3)) and the gauge invari- 
ance condition given by eq. (6). We would like to point 
out that this equilibrium spin current would exist even 
in two terminal setup. The equilibrium spin current are 
carried by all the occupied state at a given temperature. 
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This is a non-linear response and is different from linear 
response which give rise to non- equilibrium spin currents 
and is a Fermi surface property. So strictly speaking 
for the equilibrium spin currents one should take into 
account the energy dependence of spin-resolved trans- 
mission coefficient. A detailed study of the equilibrium 
spin currents would be presented in a separate article 
[16]. Here in this study we concentrate more on the non- 
equilibrium pure spin currents and the related electrical 
effects. Slonczewski has shown in Ref. [17] for magnetic 
multilayers equilibrium spin currents causes non-local ex- 
change coupling. The important difference in our case is, 
we do not need Ferromagnetic contact to have equilibrium 
spin currents, which was the case in Ref. [17]. Rather 
in our case equilibrium spin currents are generated in- 
trinsically due to the spin-orbit interaction without any 
magnetic element in the system. 

Non-equilibrium spin currents: To study non- 
equilibrium spin currents ,let us consider the case where 
the voltages at terminal 1 and 2 are respectively Vi=0 
and Vi and terminal third is a voltage probe, i.e, J| = 0. 
With this condition one can determine the voltage, Vjj, at 
third terminal using the set of equation (4) and is given 
by [10], 

Yi = _En_ (7) 

V2 T13 + T23 

The spin current flowing through terminal 3 is 
5 X £ { ( T iT - 7\V a + T 2 T - T 2 %-° )V 3 

a 

+(T^-T^)V 2 }. (8) 

From above equation (8) it is clear that 7| is non- 
zero while /| is zero by definition. Hence in terminal 3 
there is a net spin currents flowing in the absence of any 
net charge current. This is pure spin current and is in- 
trinsically generated by the spin-orbit interaction in the 
absence of any magnetization as discussed in the intro- 
duction. 

To obtain quantitative results we perform numerical 
simulation on a Y-shaped conductor shown in Fig. 1. 
We model the conductor on a square tight binding lat- 
tice with lattice spacing a and we use the corresponding 
tight binding model including spin orbit interaction given 
by eq.(l) [6]. For the calculation of spin resolved trans- 
mission coefficient , we use the recursive green function 
method. Details of this can be found in Ref. [6], [14]. 
The numerical result presented are exact and takes the 
quantum effect and multiple scattering into account. For 
the model of disorder we take Anderson model, where 
on-site energies are distributed randomly within [-U/2, 
U/2], where U is the width of distribution. All the cal- 
culation were performed on Y-shaped device of width 
d—20a, where a is lattice spacing. Other parameters are 
given in figure captions. 
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FIG. 2. spin current given by (eq.(8)) and potential given 
by (eq.(7)) versus quantization axis for different values of dis- 
order potential strength shown in inset. Calculations were per- 
formed on a device width d = 20a (see Fig.l) ,fcja = 1 , 
dimensionless spin-orbit parameter A/a 2 = 0.05 

In Fig. 2 we show the spin currents if flowing through 
the terminal 3 (Right panel) and the corresponding volt- 
age V3 when all the three terminal are non-magnetic. In 
Fig. 2 , 9=0 corresponds to z axis and 0=90 corresponds 
to y axis, we have kept fixed <fi = 90. We see that the 
maximum amount of spin currents flow along z axis. This 
is understandable since for strictly two dimensional case , 
the spin-orbit coupling given by eq. (1) conserves z com- 
ponent of spin. Hence the asymmetric scattering pro- 
duced by spin orbit interaction causes a pure spin cur- 
rents along z axis, as discussed in introduction. For the 
ballistic case (curve for U / Ep = 0) , spin currents is zero 
since there is no spin-orbit interaction in this case, as can 
be seen from eq. (1) by putting the potential U(r) = 0. 
Also for strong disorder spin current changes sign (curve 
for U/Ep = 2) due to multiple scattering. We would 
like to stress that by definition charge current flowing in 
terminal 3 is zero, i.e. /| = 0. Now from the right panel 
in Fig. 2 we see that the voltage V3 measured is differ- 
ent although there is no charge current flowing and the 
magnitude of V3 is directly proportional to the z com- 
ponent of spin current. As is seen , with the increase 
of disorder strength the magnitude of spin currents in- 
creases and accordingly the potential V3 also increases. 
However potential V3 is independent of quantization axis 
since the voltage probe is non-magnetic. Hence with a 
non-magnetic voltage probe one can detect the spin cur- 
rent , but can not measure it. To measure the spin cur- 
rents one would need a ferromagnetic voltage probe. An 
intuitive understanding of this can be gained as follows. 
From Fig. 2 (right panel) we notice that the spin currents 
depends on the quantization axis. Thus if the probe is 
a ferromagnetic, electrons which are polarized parallel 
to the ferromagnet would be transmitted easily than the 
electrons polarized anti-parallel to the ferromagnet. Since 
the voltage at the probe is determined by the ration of 
transmission coefficient (eq.(7)), hence the probe voltage 
should show variation in phase with the spin currents. 

This is confirmed in Fig. 3. Where we have plotted 
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spin current (left panel) and voltage (right panel) for the 
case when the third terminal is a ferromagnetic. Left 
panel shows the spin currents and the corresponding volt- 
age is shown in right panel. The quantization axis is 
given by the direction of magnetization. We see that as 
the spin current changes the corresponding voltage mea- 
sured also changes in phase. Hence by having a ferro- 
magnetic voltage probe one can measure the pure spin 
current. We would like to mention that in our numer- 
ical simulation voltage probe is an invasive one, i.e. it 
is strongly coupled to the system, hence one sees the 
quantitative difference between the results of Fig. 2 and 
Fig. 3. In Fig. 3 , we see that spin currents are non zero for 
the ballistic case ((U/Ep = 0)) while for non-magnetic 
case shown in Fig. 2 spin current for ballistic case is zero. 
This is so because the Ferromagnetic probe is strongly 
coupled (invasive probe), it essentially injects a polarized 
current. However this is not a hindrance for measuring 
spin currents generated by spin-orbit interaction. Since 
as is seen from Fig. 3 , it only gives rise to a constant shift 
compared to the non-magnetic case (Fig2). Recently it 
was pointed out in ref. [18] that in magnetic bilayer sys- 
tems dynamic exchange coupling arises due to the injec- 
tion of polarized charge current. In the said work effect of 
spin-orbit interaction were not taken into account. Since 
here we point out the existence of spin currents (equi- 
librium and noncquilibrium) due to the spin-orbit inter- 
action, hence such currents in principle would modify 
quantatively proposed dynamic exchange coupling. Since 
charge transport for the Y-shaped mesoscopic junction 
have been studied in past experimentally as well theoret- 
ically. In view of this we hope the study presented here 
for the spin transport would open up new opportunity in 
the field of spintronics. 
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FIG. 3. spin current flowing through the Ferromagnetic 
terminal 3 (voltage probe) versus quantization axis (left 
panel) and right panel shows the corresponding voltage at 
the terminal third. Different curves corresponds to different 
strength of disorder. Inset shows the strength of Anderson 
disorder. Ferromagnet is modeled as exchange split with ex- 
change splitting (A) given as A/Ef = 0.5. Other parameters 
are same as for Fig. 2. 



[1] S. A. Wolf,D. D. Awschalom, R. A. Buhram, J. M. 

Daughton.S. Von. et. al., Science 294, 1488 (2001) and 

references therein. 
[2] L. W. Molenkamp et al, Phys. Rev. B. 64, R121202 

(2001). 

[3] S. Datta and B. Das, Appl. Phys. Lett. 56, 665(1990). 
[4] J. E. Hirsch, Phys. Rev. Lett. 83, 1834 (1999), Shufeng 

Zhang, Phys. Rev. Lett. 85, 393 (2000). 
[5] G. Bergman, Phys. Rev. B. 63, 193101-1 (2001). 
[6] T. P. Pareek and P. Bruno, Phys. Rev. B.63, 165424-1 

(2001). P. Bruno Phys. Rev. Lett. 79, 4593 (1997). 
[7] L.D. Landau and E. M. Lifshitz Quantum Mechanics Vol. 

3, pp. 583-588. 
[8] A. A. Kiseleve and K. W. Kim cond-mat 0203261. 
[9] M. Biittiker, IBM J. Res. Develop. 32 63 (1988). 
[10] T. P. Pareek, Sandeep K. Joshi and A. M. Jayannavar 

57, 8809 (1998). 
[11] Yu. A. Bychkov and E. I. Rashba, Sov. Phys. JETP Lett. 

39, 78 (1984). 

[12] A. Brataas, Yu. V. Nazarov and G.E.W. Bauer, Phys. 

Rev. Lett 84, 2481 (2000). 
[13] M. J. Stevens, A. L. Smirl, R. D. R. Bhat, A. Najmaie, 

J. E. Sipe,and H. M. van. Driel, Phys. Rev. Lett. 90, 

136603-1 (2003). 
[14] T. P. Pareek, Phys. Rev. B. 66, 193301 (2002). T. P. 

Pareek and P. Bruno, Phys. Rev. B. 65, 241305 (2002). 
[15] S. Datta, Electronic transport in mesoscopic systems, 

(Cambridge University press, Cambridge, 1997). 
[16] T. P. Pareek, manuscript under preparation. 
[17] J. C. Slonczewski, Journal of Magnetism and Magnetic 

Materials. 126, 374 (1993). 
[18] B. Heinrich, Y. Tserkovnyak, G. Woltersdorf, A. Brataas, 

R. Urban and Gerrit E. W. Bauer, Phys. Rev. Lett. 90, 

187601-1 (2003) 



4 



